Temperature of projectile like fragments in heavy ion collisions 
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A model in which a projectile like fragment can be simply regarded as a remnant after removal of 
some part of the projectile leads to an excited fragment. This excitation energy can be calculated 
with a Hamiltonian that gives correct nuclear matter binding, compressibility and density distribu- 
tion in finite nuclei. In heavy ion collisions the model produces a dependence of excitation energy on 
impact parameter which appears to be correct but the magnitude of the excitation energy falls short. 
It is argued that dynamic effects left out in the model will increase this magnitude. The model can 
be directly extended to include dynamics but at the expense of increased computation. For many 
calculations for observables, a temperature is an easier tool to use rather than an excitation energy. 
Hence temperature dependences on impact parameter in heavy ion collisions are displayed. 

PACS numbers: 25.70Mn, 25.70Pq 



I. INTRODUCTION 

In recent times we proposed a model [IH1] for pro- 
jectile fragmentation whose predictions were compared 
with many experimental data with good success. The 
model has three parts. To start with we need an abra- 
sion cross-section. This was calculated using straight line 
trajectories for the projectile and the target leading to a 
definite mass and shape to the projectile like fragment 
(PLF). The PLF will have an excitation energy. It was 
conjectured that this will depend upon the relative size 
of the PLF with respect to the projectile, i.e., on (A s /A ) 
where A s is the size of the PLF and A$ is the size of the 
whole projectile. Since the size A s of the PLF depends 
upon the impact parameter of the collision, dependence 
on (A s /Aq) means the excitation energy depends upon 
the impact parameter 0, H. In our original version this 
dependence was neglected |l[ but became an important 
feature in the improved model 0. Instead of excita- 
tion energy we use temperature T. The hot PLF will 
disintegrate into different composites which can be cal- 
culated using a canonical thermodynamic model (CTM) 
Evaporation from hot composites which result from 
CTM was implemented 

In our model the temperature of the PLF was not cal- 
culated, it was fitted from data • I n this paper we try 
to estimate T from a more basic approach. Temperature 
T in the PLF has been studied a great deal in the past 
using a combination of theory and experimental data. A 
very popular method uses experimental populations of 
excited states (for example, the "Albergo" formula @) 
to deduce a temperature. Many data suggest that the 
temperature is of the order of 5 MeV. If the excitation 
energy in the PLF is mostly due to its crooked shape at 
the time of separation, one can estimate the excitation 
energy assuming a liquid drop model with a volume term, 
a surface tension term and Coulomb contribution. One is 
probably confined to assuming a constant density which 
is not realistic. 

The objective of this paper is to estimate the tem- 



perature in the PLF using a more microscopic approach. 
The lowest order approximation we make to estimate the 
temperature arises from the crooked shape that results 
from straight line geometry for cuts. The excitation en- 
ergy originates from structure effects solely as in a liquid 
drop model. But we are not using a liquid drop model 
nor are we limited to constant density. Our method can 
be extended to consider dynamic effects using a trans- 
port model, specifically the BUU (Boltzmann-uehling- 
Uhlcnbeck) model. The techniques we use, even for the 
lowest order estimates, are well-known in BUU calcula- 
tions. 

We do calculation for Sn on Sn, Xe on Pb, Ni on 
Be and Ta as we have used our model for these cases 
before [IH3 and experimental results are known 043- 
In heavy ion collisions we consider one ground state 
nucleus hitting another ground state nucleus. We first 
turn to the model of the ground state. 



II. MODEL OF THE GROUND STATE:- 

We use Thomas-Fermi (TF) solutions for ground 
states. Complete details of our procedure for TF solu- 
tions plus the choice of the interactions are given in Ref. 
[i~ol |. For completeness the prescription is outlined. The 
kinetic energy density is given by 



T(r) = J d 3 pf(r,p)p 2 /2m 



(1) 



where f(r,p) is the phase space density. Since we are 
looking for lowest energy we take, at each r, f(r,p) to be 
non-zero from to some maximum pp(r). Thus we will 
have 



f(r,p) 



h 3 



0\p F (r,p) -p] 



(2) 



The factor 4 is due to spin-isospin degeneracy and us- 
ing the spherical symmetry of the TF solution we have 
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dropped the vector sign on r and p. This leads to 




For potential energy we take 

V = Aj d 3 r^ft + -L.B f p° +l {r)d\ 

+- J d 3 rd 3 r'v(f,f')p(r)p(r') (4) 

The first two terms on the right hand side of the above 
equation are zero range Skyrme interactions. The third 
which is a finite range term is often suppressed and the 
constants A,B,a are chosen to fit nuclear matter equi- 
librium density, binding energy per nucleon and com- 
pressibility. In heavy ion collisions, for most purposes, 
this will be adequate but for what we seek here, possi- 
bly a small excitation energy, this is wholly inadequate. 
Thomas- Fermi solution is obtained by minimising T + V. 
With only zero range force, p(r) can be taken to be a 
constant which goes abruptly to zero at some rg fixed by 
the total number of nucleons. Now if p is chosen to min- 
imise the energy then, a nucleus, at this density with a 
cubic shape is as good as a spherical nucleus. Besides the 
minimum energy nucleus will have a sharp edge, not a re- 
alistic density distribution. This problem does not arise 
in quantum mechanical treatment with Skyrme interac- 
tion. Including a finite range potential in TF one recovers 
a more realistic density distribution for the ground state 
and one regains the nuclear structure effects which will 
contribute to excitation the PLF. This is discussed in 
more detail in Ref. [10|. 

We note in passing that Lenk and Pandharipande in- 
troduced a diffuse surface by modifying the kinetic energy 
term [14] . 

Thomas-Fermi solutions for relevant nuclei were con- 
structed with following force parameters. The constants 
A, B, and a (Eq.4) were taken to be A=-1533.6 MeV 
fm 3 , B=2805.3 MeV fm 7 / 2 , a = 7/6. For the finite range 
potential we chose an Yukawa :V y . 

e -\r-r'\/a 



with Vb=-668.65 MeV and a=0.45979 fm. Binding 
energies and density profiles for many finite nuclei with 
these parameters (and several others) are given in Ref. 
10]. These have been used in the past to construct TF 
solutions which collide in heavy ion collisions [ll| . 



III. METHODOLOGY:- 

We use the method of test particles to evaluate ex- 
citation energies of a PLF with any given shape. The 
method of test particles is well-known from use of BUU 



models for heavy ion collisions [12J. Earlier applications 
were made by Wong [13[ . 

We first construct a TF solution using iterative tech- 
niques [13 . The TF phase space distribution will then 
be modeled by choosing test particles with appropriate 
positions and momenta using Monte-Carlo. Throughout 
this work we consider 100 test particles (N tes t = 100) 
for each nucleon. For example, the phase space distribu- 
tion of 58 Ni is described by 5800 test particles. A PLF 
can be constructed by removing a set of test particles. 
Which test particles will be removed depends upon col- 
lision geometry envisaged. For example, consider central 
collision of 58 Ni on 9 Be. Let z to be the beam direction. 
For impact parameter 6=0 we remove all test particles in 
58 Ni whose distance from the center of mass of 58 Ni has 
x 2 + y 2 < rg 2 where rg = 2.38 fm is the radius at half 
density of 9 Be. The cases of non-zero impact parameter 
can be similarly considered. 

The " sudden approximation" we consider is the follow- 
ing. We assume that the PLF is formed suddenly. At the 
time the PLF separates from the participants the shape 
and momentum distribution of the PLF can be described 
by removing some test particles as described above. Of 
course this PLF will undergo many more changes later 
but all we are concerned with is the energy of the system 
at the time of "separation". Since the PLF now is an 
isolated system, the energy will be conserved. Of course 
the Coulomb force from the participants will continue to 
be felt by the PLF. But the major effect of this will be on 
overall translation of the PLF and all we are interested 
in is intrinsic energy. 

We now describe how we calculate the energy of this 
"crooked" shape object. The mass number of the PLF 
is the sum of the number of test particles remaining di- 
vided by Ntest- Similarly the total kinetic energy of the 
PLF is the sum of kinetic energies of the teat particles 
divided by N tes t- Evaluating potential energy requires 
much more work. We need a smooth density to be gener- 
ated by positions of test particles. We use the method of 
Lenk and Pandharipande to obtain this smooth density. 
Other methods are possible [l2|. Experience has shown 
that Vlasov propagation with Lenk-Pandharipande pre- 
scription conserves energy and momenta very well [l4j 
although we will not get into time propagation in this 
paper. 

The configuration space is divided into cubic lattices. 
The lattice points are I fm apart. Thus the configuration 
space is discretized into boxes of size ^ 3 fm 3 . Density at 
lattice point r a is defined by 

AN tcBt 

Pi,(r a ) = S (r<x-ri) (6) 

i=l 

The form factor is 

S W = M l (ril \^ )g{y)g{z) (7) 
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where 

g(q) = (nl-\q\)Q(nl-\q\) (8) 

The advantage of this form factor is detailed in [l4| so 
we will not enter that discussion here. In this work we 
have always used l—l fm and n=l. 

It remains to state how we evaluate the potential en- 
ergy term (eq.(4)). The zero range Skyrme interaction 
contributions are very simple. For example the first term 
is calculated by using 

A [ d 3 r^P- = Aj2(l 3 )pl(r a )/2 (9) 

J a 

With our choice Z 3 =l fm 3 . For the third term in eq.(4) 
(the Yukawa term) is rewritten as 1/2 ^2 a PL('f r a)4'L('r'~a) 
where 4>{r) is the potential at f due to the Yukawa, 
i.e., 0(f) = JVy(\r-P\)p(r')d 3 r'. 

The calculation of Yukawa (and/or Coulomb) potential 
due to a charge distribution which is specified at points 
of cubic lattices is very non-trivial and involves iterative 
procedure. This has been used a great deal in applica- 
tions involving time-dependent Hartree-Fock theory. We 
will just give the references fl5l - [T7j |. We also found an 
unpublished MSU report very helpful [llj]. 

With this method we can calculate the total energy of 
the PLF. However we are interested in excitation energy 
of the system which requires us to find the ground state 
energy of the PLF which has lost some nucleons. We can 
use TF theory to find this. The iterative TF solution 
also gives the ground state energy. But then we will be 
using two different methods for evaluating energies, one 
for the PLF crooked shape and a different one for the PLF 
ground state. Although the results are quite close, it is 
more consistent to use the same prescription for ground 
state energy and excitation energy. Hence, even for the 
ground state energy we generate test particles and go 
through the same procedure as for the PLF with crooked 
shape. 

However we are not finished yet. This gives us the 
excitation energy but we need to know the temperature 
corresponding to this excitation. The canonical thermo- 
dynamic model (CTM) [J] can be used to calculate aver- 
age excitation per nucleon for a given temperature, mass 
number and charge. This model is used to deduce the 
temperature for the PLF from its excitation energy. Re- 
sults are shown in the next section. 



IV. RESULTS:- 

We show results for 124 Sn on 119 Sn, 136 Xc on 208 Pb, 
58 Ni on 9 Be and 181 Ta as these were studied in our pre- 
vious work. We used the data of Sn on Sn in conjunction 
with our model to deduce " experimental" values of tem- 
perature as a function of measured PLF mass . For Sn 
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FIG. 1: Plot of temperature against mass number of the PLF 
for 124 Sn on 119 Sn reaction. The solid curve is the present 
model, the dotted line is from our previous work 0]. Since 
in this particular case the previous work fixed the tempera- 
ture by matching with data on intermediate mass fragment 
production, the dashed curve represents, in some sense, "ex- 
perimental" data. The dashed curve was parametrised as a 
function of the wound of the PLF defined as 1.0-A 3 /Ao- The 
same parametrisation was used for dashed curves in Figures 
3, 5 and 6. 



on Sn, the total measured charge in the forward direction 
specified the size of the PLF. In the calculation we can 
match this size by varying the impact parameter. For 
this size, experiment measured the number of intermedi- 
ate mass fragments (defined to be composites with charge 
between 3 and 20). In our model, for a given size, the 
number of intermediate mass fragments depends on the 
temperature of the PLF. We exploit this to deduce the 
temperature for this b. Thus we are able to , for Sn on 
Sn, a temperature vs b curve. However a much more use- 
ful parametrisation is in terms of the wound (=1-A s /Aq) 
as this can be used for other colliding pairs. We assumed 
a universal curve derived from the experimental curve of 




Impact parameter (fm) PLF Mass (Ag) 

FIG. 2: In the model of this paper, plot of excitation energy 
per nucleon for 124 Sn on 119 Sn against (a) impact parameter 
and (b) mass A s of the PLF. 
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Sn on Sn and used that to calculate many cross-sections 
(not only intermediate mass fragments) in all the cases 
mentioned above. The agreements with data were very 
pleasing. In our present case, no adjustable parameters 
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FIG. 3: Curve similar to that of fig.l but for 136 Xe on 208 Pb. 
The dotted curve is from our previous work. 

are available to deduce either the excitation energy or the 
temperature. All we have is a semi-classical Hamiltonian 
which reproduces the ground state energy, correct radius 
at half density and a realistic distribution of density in 
finite nuclei. It fits nuclear matter density, binding en- 
ergy and compressibility. Consider the projectile 124 Sn 
hitting the target 119 Sn. We obtain the TF ground state 
of the projectile and create test particles with appropri- 
ate positions and momenta which map this phase-space 
distribution. The effect of collision with the target is that 
some test particles of the projectile are driven out leaving 
a PLF with smaller mass number and a crooked shape. 
A disc with the radius of the target will take out different 
numbers of projectile test particles for different impact 
parameters. We have already described how excitation 
energies are calculated. 

Figures 1 and 2 show the case of 124 Sn on 119 Sn. In Fig- 
ure 1 the curve showing previous calculation is essentially 
the experimental result as the temperatures which repro- 
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FIG. 5: Same as in Fig. 3 except for the case of 5S Ni on 9 Be. 



duce experimental results were selected. Our present re- 
sults obtained solely from a Hamiltonian which fits just 
ground state data does have the correct trend but un- 
derestimates the temperature. This is to be expected. 
We have assumed zero coupling between participants and 
spectators. Participants have much more energy density 
(hence higher temperature). Energy flows from higher 
value to lower value and hence temperature of the spec- 
tator is expected to rise further. Figure 2 plots excitation 
energy per nucleon against impact parameter and mass 
number of the PLF. In this paper excitation energy is 
directly calculated. Temperature is deduced from excita- 
tion energy using the CTM. 
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FIG. 4: Same as in Fig.2 except for 136 Xe on 208 Pb. 
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FIG. 6: Same as in Fig.5 except for 58 Ni on 181 Ta. 

Figures 3 and 4 show similar results for Xe on Pb. 
The functional dependence of the temperature on the 
"wound" is taken from Sn on Sn. Previous calculations 
gave good fits with many experimental observables but 
they are not displayed here and can be found in ref. [2, 3] . 
Figures 5 and 6 show our results for 58 Ni on 9 Be and 
on 181 Ta. The results for Ni on Ta are similar to those 
of Sn on Sn and Xe on Pb. For Ni on Be the present 
model makes the fall of temperature with growing mass 
A s much faster. This may have something to do with the 
rather small size of the target. For the smallest possible 
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target, i.e., a proton, the model we have been following 
here is clearly inappropriate. At what size of the target, 
the model becomes reasonable is hard to estimate. 

The model and the methodology used here can be 
directly extended to BUU calculations to examine the 
approximations and possible improvements. We have 
embarked on such a project. Preliminary estimates are 
promising. 
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